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Residual stress are stresses that remains in an engineering material, 
especially metallic materials, after the original applied load has been 
removed. Its measurement will help reveal areas in fabricated materials 
highly stressed. This paper presents a design of a portable residual stress 
measuring device based on the magnetic barkhausen noise principle. This 
was achieved by using two permanent magnets, a pickup coil, amplifier 
and an analogue to digital converter to create the hardware. The software 
was designed using LabVlEW 2019. After fabrication, measuring 
experiments on mild steel materials were conducted and the results 
obtained from using the device for measurement were compared to 
simulated values (Simufact.Weld 6.0) CAD software. Comparing its 
results with that obtained from Simulation showed an R^ value of 77.47% 
with a strong correlation of 0.83land a P-value of0.0071 which was less 
than the 0.05 alpha level taken. This device can be adopted in weld 
residual stress measurement by artisans and roadside welders for 
effective minimization of residual stress and failures associated with it. 


1. Introduction 

Nearly all infrastructural disasters recorded in Nigeria are manmade [1]. Technical and mechanical 
failures of constructed and fabricated materials which pose danger to live and property and, in 
some cases, claimed lives is on the increase. These failures are seen in a vast range of products 
including collapsed buildings, broken bill board stands, water storage stand failure, toppling trailer, 
refurbished automobile part failure, even house hold appliances [2]. Most fabrication works in 
Nigeria are carried out by road-side artisans who are ignorant of the danger posed by residual 
stress. Residual stress are stresses that remains in an engineering material, especially metallic 
materials, after the original applied load has been removed [3] and [4]. These stresses are usually 
introduced into the materials, during manufacturing, fabrication or during machine operation [5]. 
Residual stresses are measured using destructive, semi-destructive or non-destructive methods [6], 
the finite element method (FEM) are seldom used for residual stress measurement because it is 
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time consuming. Manufactured and fabricated components are best measured using the non¬ 
destructive method as finished products will not be destroyed in the process. X-ray diffraction 
technique is the most popular method of non-destructive residual stress measurement, and 
measurement is often carried out in the lab one major drawback is that the specimen is limited by 
the size of the device x-ray chamber. In Nigeria, these measuring devices are hard to find and very 
expensive to purchase by researchers and fabricator artisans. Most of the structural failures 
experienced in Nigeria, can be minimized or prevented if an affordable residual stress safety device 
can introduce to the Nigerian fabrication industries, since more ferromagnetic (mostly iron and its 
alloys) materials are mostly used for welding. Residual stress measurement will help reveal areas 
in fabricated materials highly stressed, it also aid designers in estimating good shelf life and also 
minimizes “unexpected” failures and the unpleasant situations associated with them. 

This study presents a design of a portable residual stress measuring device, based on the magnetic 
barkhausen noise, that can be used as an alternative to X-ray diffraction method for ferromagnetic 
materials and also to demonstrate its application in welding as a promising weld residual stress 
measuring device which can produce a fast measurement indoor and outdoor. 

2. Materials and Methods 

This device was designed based on the magnetic barkhausen noise principle. Permanent magnets, 
search coil (also known as pickup coil), mini amplifier (amp), analogue to digital converter (ADC) 
and a software to analyze the data (Labview), were the materials used for design. 

Two identical strong rare earth neodymium iron boron (NdFeB) magnet (N52) were used for the 
magnetic sensor. The magnet is a cylinder with diameter 10mm and height of 20 mm. These 
magnets are brittle, therefore, extreme care was taken to avoid damaging it by not allowing it 
smash together due to its attractive power. 

Equation (1) given by [7] help us understand the formation of magnetic domains as it occurs 
because of a minimization contest of the five basic energies involved in ferromagnetism: 

E — fi'exc/ian^e Emagnetostatic Emagnetocrystalline Emagnetoelastic E^all (1) 


The penetration depth of magnet on ferromagnetic material can be calculated (estimate) using the 
formula obtained from [8], presented in Equation (2). 

6 = 1/V (nuaf) (2) 

Where 5 represents the penetration depth, p denotes the magnetic permeability, a represents the 
electrical conductivity while / represents the frequency of the alternating magnetic field. 

In other to determine the depth of magnetic influence on our material, the experimental (crude 
method) means was adopted. Firstly, we started by collecting some iron dust from Engineering 
Laboratory (ELA), University of Benin. We considered six (6) categories of specimen thicknesses 
for the experiment namely: 1mm plate, 2mm plate, 3mm plate, 4mm plate, 5mm plate and 6mm 
plate which are all ferromagnetic mild steel materials. 
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Figure 1: Illustration of the depth of penetration experiment 


Air gap of 2mm between the magnet and the plates was created, with the magnet placed on the 
plate while the plate was placed of the iron dust as sketched in Figure 1. Table 1 gives the summary 
results observed from the experiment. 


Table 1. Depth of penetration 


Mild steel Plate 

Magnetic influence 

observed 

N 

VL 

L 

S 

VS 

1mm 






2mm 



V 



3mm 






4mm 

V 





5 mm 

V 





6mm 

V 






Where N represent no visible effect on iron dust magnetization, VL denotes very low 
magnetization effect on iron dust, L means low magnetization effect on iron dust, S represent 
strong magnetization effect on iron dust and VS means very strong magnetization effect on iron 
dust. 

The search coil is a solenoid fitted with two permanent neodymium magnets, the induction coil 
sensor has a transfer function of F = /(S) resulting from the Faraday’s law of induction where 

V =-N — = -NA— = -uNA — (3) 

dt dt ^ dt ^ ^ 

The aim of the coil design is to obtain a strong magnetic signal using the faradays law of induction, 
which from proof of principle, the output signal, V, of a coil sensor depends on the rate of change 
of flux density, dB/dt. It was observed that the induced voltage V from an alternating magnetic 
field, is affected by the core material where the permeability (/r) is presented in equation 3. 

In the coil design, coil thickness of 1mm was used due to limited technology in handling the 
previous 0.01, 0.02 and 0.03mm coil earlier purchased for the project, as the 0.01-0.03mm coil 
wire kept breaking during winding. [10] Suggested smaller wire diameter for better results, a 
greater number of turns of wire in the coil also contributes to greater inductance; while a fewer 
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numbers of turns results in less inductance. The pickup coil had a winding of 2500 turns which 
was achieved by using a locally made winder. Figure 2 shows the completed pickup coil. 



Figure 2: pickup coil with two N52 magnets fixed in place 


The signal generated due to domain flipping in the ferromagnetic material picked-up by the search 
coil came out very faint, therefore, it was boosted/magnified with the PAM8403. PAM8403, 3W, 
class-D audio amplifier. 

To analyze the signal generated, the Arduino Nano with Atmega328p processor was used. 

The resolution provided by the Arduino Nano ADC is 10 bits, the ADSl 115 which provides a 16- 
bit resolution and communicates via I2C interface was used to enhance the sensor for better 
resolution. In summary, the ADSl 115 was used as the primary analog to digital converter (slave), 
connected to the Arduino Nano (Master). 


In the design presented in Figure 3, two of N52 magnets where employed with their dipole 
vertically opposite. The reason for that is aimed at eliminating trapping/pinning of domain walls 
that results whenever the magnet glides across the surface of the specimen, it is expected that the 
domain wall should flip due to the diploe directions which help reveal more details about the nature 
of the internal material. 
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Figure 3: Schematic image of barkhausen sensor, measuring specimen by moving across 

specimen surface with a velocity V 
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The two wires from the coil were connected to the PAM8403 mini amp, the output of PAM8403 
was connected to the ADSl 115 ADC. The ADC was connected to the nano, then to the computer. 
To establish communication between the nano and the ADC, the following connections where 
made. Note Vdd = positive voltage input to ADS 1115, GRD = ground or the negative voltage 
connection, SCL = serial clock line, SDA = serial data address 

Vdd from the ADC was connected to the 5v pin on the Nano. 

GRD from the ADC was connected to the GRD on the Nano. 

SCL from the ADC was connected to A5 on the Nano. 

Lastly, SDA from the ADC was connected to A4 on the Nano. 

For compatibility, the pickup coil, amp and ADC were fitted in a small rectangular box designed 
to be handy and have four mini tyres at the bottom to enable dynamic motion across the surface of 
the sample. 

Orthographic view for our sensor casing is shown in Figure 4, HP laptop power pack of 
50x30x125mm presented in Figure 5 was recycled, its internal components discarded and the 
casing used to house our magnets, coil, ADSl 115 and Arduino board. The four (4) tyres fitted to 
the device, had an outer diameter of 1mm. An area of 880mm^ was created at the base of the casing. 
This area serve as a vacuum through which the sensor can interact with the specimen during 
measurement, picture of the fabricated device with its connectivity is presented in Figure 6 


Top view 




T 


125m 
Left view 


1 





Right view 


Front view 


Figure 4: Orthographic view for our sensor casing 
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Figure 5: Recycled laptop power pack fitted with tyres 



Figure 6: device measurement connection 
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Lab VIEW 2019 was used to analyze the data obtained. The user interface in Figure 7 was designed 
by using sets of tools and objects. After building the front panel, we added code using graphical 
representations of functions to control the front panel objects. The graphical code, also known as 
G code or block diagram code was added, to the block diagram which resembles a flowchart. The 
block diagram, front panel, and graphical representations of code compose a VI. 


Control Measurement Sample Measurement 


Raw Volts jo 
Maximum 2 119,6625 
Minimum 2 [o 
Arithmetic Mean 2 j 8.36891 
Variance j 46.3319 
Standard Dev J 6.80675 
Skewness J.q.IOI ^ 
Kurtosis 11.56673 
signal pt per m.s |l024 


IS Volt 110.26 


RMS 10.7854 


129.8933 

|12.68 



Fault/result summary panel 

alternative hypothesis 
^mean(pop) 


50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 



confidence interval 


18.074039 

jlnf 

;test info 


sample me 

an sample std 

18.36891 

16.80675 

sample z value z value 

[1.64485 

1-24.0669 


1.000000 

reject null hypothesis 

m ^ 

W 

signifKance level 


I std SignifKance leve 

- I"” ir 


Figure 7: Lab VIEW GUI for the Magnetic device 


To validate our result, comparison between results obtained between repeated experimental weld 
reading for consistency and simulated values obtained using simufact.weld 6.0 software in Figure 
8 were conducted, using weld parameter in Table 2 . 


Table 2: Weld Process Parameter 


Sample 

Current (Amp) 

Voltage (V) 

Weld Speed (mm/s) 

Bll 

120 

62 

2.497 



Figure 8: Yield stresses for specimen Bu before weld 


76 





















































C.E Etin-osaet al. / NIPES Journal of Science and Technology Research 
2(3) 2020 pp. 70-84 

Thirty seven points were selected from the simulated specimen as shown in Figure 8. The 
yield stresses across the material before weld is given as 307.3 N/mm^. 


3. Results and Discussion 

3.1. Comparison between device results and Simufact.weld 6.0 

In this study, five (5) repeated readings were obtained for Bu using the magnetic device. The wave 
form chart readings for Bn was exported to Excel and Minitab for analysis. One out of the five 
readings was compared with result obtained from the simulated welded component in 
Simufat.weld. 

The results obtained from the above experiment for five (5) repeated readings and from 
simufact.weld are presented in Table 3, for specimen Bu. Thirty seven (37) signals were obtained 
from the specimen for one experimental run, the five repeated measurements were compared 
amongst themselves for any significant variability. The signals for each of the repeated experiment 
against time are presented in Figure 9-13. Figure 14 shows the combine series time plot for the 
five (5) of specimen Bu, while Table 3 present the Magnetic Device reading for specimen Bu. 



Time 


Figure 9: F' reading from Bll 
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Figure 11: 3rd reading from B11 



Time 


Figure 12: 4^^ reading from B1 



Time 


Figure 13: 5* reading from B11 


Table 3: Signal results generated from specimen Bu using the fabricated magnetic device in 

millivolts 


Magnetic Device reading for specimen Bn 

S/N 

1st 

reading 

mV 

2nd 

reading 

mV 

3rd 

reading 

mV 

4th 

reading 

mV 

5th 

reading 

mV 

Simufact 

reading 

N/mm^ 

1 

12.95 

13.1501 

12.6489 

12.849 

13.0501 

0 

2 

15.1971 

15.095 

15.0261 

14.894 

14.8255 

0 

3 

16.9532 

16.6531 

16.6521 

17.0624 

17.0533 

0.1 

4 

16.7997 

16.9998 

16.9309 

17.131 

16.7303 

0.1 

5 

15.7911 

15.689 

15.49 

15.488 

15.8912 

0.2 

6 

15.5948 

15.2947 

15.2258 

15.4259 

15.0252 

0.4 

7 

15.0768 

15.2769 

14.7757 

15.0759 

15.1769 

0.7 

8 

15.8554 

15.7533 

15.6844 

15.9646 

15.9555 

1.7 

9 

16.7733 

16.4732 

16.4722 

16.6723 

16.2716 

3.1 

10 

16.9787 

17.1788 

17.1099 

16.9778 

17.0788 

6 

11 

18.0628 

17.9607 

17.7617 

17.9618 

17.5611 

10.7 
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12 

17.1557 

16.8556 

16.7867 

16.6546 

17.2558 

18.2 

13 

18.9685 

19.1686 

18.6674 

19.0777 

18.4668 

31.4 

14 

18.0254 

17.9233 

17.8544 

18.0545 

18.1255 

54 

15 

17.7502 

17.4501 

17.4491 

17.2491 

17.8503 

74.1 

16 

18.3294 

18.5295 

18.4606 

18.6607 

18.26 

94.1 

17 

19.3058 

19.2037 

19.0047 

19.0027 

19.4059 

121.3 

18 

19.1585 

18.8584 

18.7895 

19.2677 

18.5889 

129.2 

19 

19.0515 

19.2516 

18.7504 

18.9505 

19.1516 

149.1 

20 

19.358 

19.2559 

19.187 

19.0549 

18.9864 

140.9 

21 

18.3155 

19.0154 

18.0144 

18.2145 

18.4156 

127.7 

22 

18.4627 

18.6628 

18.5939 

18.4618 

18.5628 

109.9 

23 

17.5904 

17.4883 

17.2893 

17.6996 

17.0887 

83.5 

24 

16.5015 

16.2014 

16.1325 

16.3326 

16.6016 

64.9 

25 

16.4052 

16.6053 

16.1041 

16.4043 

15.9035 

42.1 

26 

17.2461 

17.144 

17.0751 

17.2752 

17.3462 

25.3 

27 

16.177 

15.8769 

15.8759 

15.6759 

15.6753 

14.5 

28 

16.2796 

16.4797 

16.4108 

16.3888 

16.3797 

8.5 

29 

17.3237 

17.2216 

17.0226 

17.2227 

17.4238 

4.4 

30 

15.3248 

15.0247 

14.9558 

14.8237 

14.7552 

2.3 

31 

15.5936 

15.7937 

15.2925 

15.4926 

15.6937 

1.2 

32 

15.7148 

15.6127 

15.5438 

15.4117 

15.3432 

0.6 

33 

15.3768 

15.0767 

15.0757 

15.486 

15.4769 

0.3 

34 

15.4212 

15.6213 

15.5524 

15.7525 

15.3518 

0.2 

35 

15.2856 

15.1835 

14.9845 

14.9825 

15.3857 

0.1 

36 

15.1853 

14.8852 

14.8163 

15.0164 

15.2854 

0 

37 

14.4525 

14.6526 

14.1514 

14.4516 

13.9508 

0 


Time Series Plot of 1st reading, 2nd reading, 3rd reading, ... 



Variable- 
1st reading 
2nd reading 
Jrd reading 
4-tb reading 
Stb reading 


Figure 14: Time series plot of combined reading for Bu 
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The yield stresses on the surface of the simulated material in Figure 15 is given a default colour of 
dark red across the surface of the material before weld. During weld, this material yield stresses of 
307.3 N/mm^ becomes affected due to heat input from the weld source. Different colours 
symbolizes the extent of yield on the material as it undergoes plastic deformation due to heat input, 
the dark blue region presented in Figure 16 shows the area mostly affected during the welding 
process. 



Figure 15: Yield stresses before welding of specimen Bu 


The yield point is the point on a stress-strain curve that indicates the limit of elastic behavior and 
the beginning plastic behavior. Prior to the yield point, the material will deform elastically and will 
return to its original shape when the applied stress is removed. Residual stresses arise whenever a 
component is stressed beyond its elastic limit and plastic deformation occurs. Plastic deformation 
occurs when the stress exceeds a metal's yield strength. In this research, the difference between 
yield stress before and after weld was taken as the amount of residual stress in the material. This 
was achieved by subtracting the yield stress after weld from the yield stress before weld. The area 
with yield stress of 158.2 N/mm^ implies that it has accumulated residual stress of 149.1 N/mm^. 



Fig. 16: Yield stresses after welding of specimen Bn 
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The residual stress developed due to the weld activity is presented in Figure 17 



4 8 12 16 20 24 28 32 36 

Index 


Figure 17: Residual stress plot for specimen Bu 

Correlation is a statistical tool that is used to show how strongly pairs of variables are related. 
We applied this statistical tool to measure how strong the relationships are between the 1st 
up to the 5th experiment and Simufact.weld. Positive correlations were observed from the 
scatter plot presented in Figure 18, which shows a good agreement of experimental and 
simulated. 


tress IM/mm2, B11 (Initl) - B11 ClFinial>lsJ/mm2, 1st readinig. 2nd reading, 3rd reading 



Figure 18: Correlation matrix for Bu experimental and simulated 

For this analysis. Table 4 shows the correlation summary of five (5) signals and the 
simufact reading. Strong positive relationship was observed for all signals generated. The 
lowest correlation of 0.831 was observed between F' and the Simufact.weld which implies 
a strong correlation. 

A quadratic model was used for our fitted line plot as presented in Figure 19 which contains 
curvature. 
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For P-value greater than alpha level, the association is regarded as not statistically significant. 
The equation for Regression Analysis for the F' reading (y) versus Simufact.weld (x) result 
is presented in equation 4. The results presented in table 6 and 7 shows the p-value for Bn to 
be 0.0071, which is less than the significance level of 0.05. These results indicate that the 
association between the experimental and simulated (Simufact.welding) is statistically 
significant. 


Table 4: Correlations Matrix for Bn Experimental and Simufact.we' 



SMF.W.Bll 

1st 

reading 

2nd 

reading 

3rd 

reading 

4th 

reading 

1st reading 

0.831 





2nd reading 

0.835 

0.993 




3rd reading 

0.862 

0.982 

0.983 



4th reading 

0.853 

0.979 

0.982 

0.995 


5 th reading 

0.861 

0.983 

0.979 

0.980 

0.978 


d 


y = 7.774 + 0.1792X - O.OOOSObx^ (4) 

Table 5 shows the model summary of 2.23562 for s, 77.47% for coefficient of 
determination and 76.14% for the adjusted R squared where S represents the average 
distance that the observed values from the regression line. Most importantly, it shows 
how wrong the regression model is on average using the units of the response variable. 
Smaller values are preferable because it indicates that the observations are closer to the 
fitted line. R^ is regarded as percentage of variation of the response that can be 
explained by the model. The higher the R^ value, the better the model fits the data. R^ 
is always between 0% and 100%. We used adjusted R^ when we compare models that 
have different numbers of predictors (Simulated vs experimental). R^ always increases 
when predictors are added to the model, even when it will not improvement to the 
model. The adjusted R^ value incorporates the number of predictors in the model 
analysis to choose a better model. 


Table 5: 

: Mode] 

Summary 

s 

R-sq 

R-sq(ad.j) 

2.23562 

77.47% 

76.14% 


Table 6: Analysis of Variance 


Source 

DF 

SS 

MS 

F 

P 

Regression 

2 

584.264 

292.132 

58.45 

0.0071 

Error 

34 

169.932 

4.998 



Total 

36 

754.196 
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Tab 


e 7: Sequential Analysis of Varia nce 


Source 

DF 

SS 

F 

P 

Linear 

1 

520.844 

78.12 

0.007 

Quadratic 

1 

63.420 

12.69 

0.001 


Fitted Line Plot 

1st reading = 7.774 + 0.1792 SMF B11 
- 0,000306 SMF 611^2 



s 

223562 

R-Sq 

77.5% 

R-Sq|adj) 

76.1% 


Figure 19: Fitted line plot of Bn simulated 


4.1. Discussion 

The use of magnet to measure residual stress signals has been sueeessfully demonstrated. The 
Pearson Correlation type was selected for this analysis, it measures the strength and direction of a 
linear relationship between two variables. The values range from -1 (strong negative relationship) 
to +1 (strong positive relationship). Values at zero or close to zero implies no or weak 
relationship/correlation. The lowest correlation between simulated reading vs experimental 
reading for Bn was found to be 0.831 which indicate a very strong correlation. 

Regression analysis was used to determine if the association between the response and each term 
in the model is statistically significant, and was also applied to determine if the 1st experiments vs 
Simufact.weld is statistically significant. The p-value obtained was compare to our significance 
level to assess the null hypothesis. The null hypothesis is that the term's coefficient is equal to zero, 
which indicates that there is no association between the term and the response. A significance level 
of 0.05 was selected. A significance level of 0.05 indicates a 5% risk of concluding that an 
association exists when there is no actual association. The P-value less than the alpha level shows 
the association to be statistically significant, and we can conclude that there is a statistically 
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significant association between the response variable and the term. For regression analysis, S is 
the average distance from the regression line. It was observed that the regression model developed 
for Bii simulation vs experimental had S values of 2.23562. Coefficient of determination (R^) 
values of 77.47% for Bn. Adjusted values of 76.14% for Bn for each simulated reading vs 
experiments model. Comparing the readings from the device to the simulated value, it was 
observed that there was a good agreement between the experimental and the simulated result. 

5. Conclusion 

A promising simple and effective residual stress measuring option, locally fabricated has been 
successfully demonstrated. Comparing its results with that obtained from simulated means, we 
obtained R^ of 77.47% with a strong correlation of 0.831. It had a P-value of 0.0071 which was 
less than the 0.05 alpha level taken. This device can be adopted in weld residual stress 
measurement by artisan and roadside welders for effective minimization of residual stress and 
failures associated with it. 
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